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Abstract 



We liave considered short-baseline neutrino oscillations, /3-decay and 
(/3/3)o!/-decay in two schemes with an inverted mass spectrum and mixing 
of three and four massive neutrino fields. We have analyzed the results of the 
latest experiments on the search for oscillations of terrestrial neutrinos and we 
have discussed the compatibility of the LSND indication in favor of neutrino 
oscillations with the results of the other experiments. In the framework of the 
models under consideration, it is shown that the observation of (/3/5)oiy-decay 



1 



could allow to obtain information about the CP violation in the lepton sector. 
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I. INTRODUCTION 



The problem of neutrino mass and mixing is the central issue of today's neutrino physics. 
The study of this problem is associated with the search of new scales in physics. The 
investigation of the neutrino mass problem is also connected with the hopes to understand 
the nature of the dark matter in the universe. 

At present there are several indications in favor of non-zero neutrino mass and mixing. 
One indication comes from the solar neutrino experiments. In all four solar neutrino ex- 
periments (Homestake Kamiokande 0, GALLEX [|| and SAGE [Q) the observed event 
rates are significantly lower than the event rates predicted by the Standard Solar Models (see 
Refs. After the test of the GALLEX and SAGE detectors in special experiments with 

radioactive ^^Cr sources 1^,^] the indications in favor of neutrino oscillations coming from 
solar neutrino experiments have become more significant. Furthermore, a phenomenological 
analysis of the solar neutrino data based on the assumption that there is no neutrino mixing 
indicates that the signals due to the ^Be neutrinos in the Cl-Ar and Ga-Ge experiments 
are substantially lower than the signals predicted by the Standard Solar Models. No plau- 
sible astrophysical or nuclear physics explanation of this discrepancy has been proposed so 
far. 

The solar neutrino data can be explained by assuming neutrino mixing with the MSW 



effect [[lT| , p!2| or vacuum oscillations [p!^-[T5||. In the simplest two-generation scheme, with 
transitions between z/g and z/^ or z/g and v^, a fit of the data with the MSW mechanism yields 
the following values for the mixing parameters |T^: Am^ ~ 5 x 10~^eV^ and sin^26' ~ 
7 X 10~^ or Am^ ~ 2 x 10~^ eV^ and sin^ 26 ~ 0.8, where Am^ = — m\ (mi and m2 are 
the neutrino masses) and 9 is the mixing angle. In the case of vacuum oscillations, a fit of 
the data yields ^ Am'^ 6 x 10"^^ eV^ and sin^ 20 ~ 0.9. 

Another indication in favor of neutrino mixing comes from the results of the experiments 
on the detection of atmospheric neutrinos [|T^ -pO|. In the Kamiokande, 1MB and Soudan 



experiments the ratio of //-like and e-like events is less than expected. The data can be 
explained by ^ Ur or z/^ ^ z/g oscillations with Am^ ~ 10~^ eV^ and a large mixing 
angle. 

A third indication in favor of neutrino mixing comes from the possible evidence of z/^ — > 
transitions reported recently by the LSND collaboration |^ (see, however, also Ref. [0). 
The explanation of the LSND data in terms of z/^ ^ z/g oscillations requires a value of Am^ 
of the order of a few eV^. 

Finally, neutrinos are very plausible candidates for the dark matter in the universe if 
their masses lie in the eV region (see, for example, Refs. p3| , ^ ). 

All these indications in favor of neutrino mass and mixing could imply that there are at 
least three different scales for Am^. In this case, the number of massive neutrinos must be 
more than three. Let us notice that the results of the LEP experiments which proved that 



the number of neutrino flavors is equal to three (see Ref. [^) do not constrain the number of 



massive light neutrinos. In the case of a Dirac neutrino mass term, the total lepton number 
is conserved, massive neutrinos are Dirac particles and the number of massive neutrinos can 
be equal to three. In the case of a Dirac and Majorana neutrino mass term, which is typical 
for GUT models, the total lepton number is not conserved, massive neutrinos are Majorana 
particles and their number is larger than three. 



3 



We will consider here short-baseline oscillations of the terrestrial neutrinos in two possible 
schemes, one with mixing of three massive neutrino fields and the other with mixing of four 
massive neutrino fields. In the case of three neutrinos we will assume that the neutrino 
masses satisfy the relation mi -C m2 — m^. In the case of four neutrinos we will assume that 
mi ~ 1712 <^ ^3 — ^i- Models with three and four massive neutrinos with this types of mass 
relations have been considered in the recent articles p6|-[28[] and [^9| -|3^ respectively. They 



are inspired by the experimental indications in favor of neutrino mixing and by astrophysical 
arguments |22[ in favor of the existence of two practically degenerate neutrinos with masses 
in the eV range. These schemes are compatible with the constraints that follow from the 



r-process production of heavy elements in the neutrino-heated ejecta of supernovae |^ . 

We will derive now the general formulas for neutrino oscillations in the framework of the 
schemes under consideration, that will be used in the subsequent analyses. Let us consider 
neutrino oscillations in the case of two groups of neutrinos with close masses. We will assume 
that the masses of the second group are much larger than the masses of the first group and 
only one squared-mass difference is relevant for oscillations of the terrestrial neutrinos. 

According to a general theory of neutrino mixing (see, for example, Refs. |I^p6| , p4| ), the 



left-handed flavor neutrino fields Vai are superpositions of the left-handed components of 
(Dirac or Majorana) massive neutrino fields z/ji 



"^Uaii^iL , a = e,fj,,T. (1.1) 



If the number of massive neutrinos is more than three, we also have 

(J^sr)" = ^UsiUiL , S = Si,S2,..., (1.2) 

i 

where {usrY = Cv^r {C is the matrix of charge-conjugation), are right-handed (sterile) 
fields, [/ is a unitary mixing matrix and the index s denotes the last (one or more, depending 
on the number of the sterile fields) rows of U. From Eqs.( |l.l|) and ( |1.2|) it follows that not 
only transitions among active neutrino states are possible, but also transitions among active 
and sterile states. We will consider the case of two neutrino groups with close masses, 

nil < ... < fnM and mM+i < • • • < rriN , (1.3) 

and we will assume that 

m,«m,, with [i = M + l,...,N. ^^'^^ 
We will also assume that in the experiments with terrestrial neutrinos we have 

(m^ — m^,) L , 



2p 

(mi — ml,) L 

^ ^ < 1 , k,k'>M , 

2p , , , 



(1.5) 
(1.6) 
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where L is the distance between the source and the detector and p is the neutrino momentum. 
Taking into account Eqs. (|1.5|) and (|1.6|) , the amphtude of — > z^/? transitions (here and 
i>j3 are the active or sterile neutrinos) for the short-basehne terrestrial experiments is given 
by 



N 



A, 



- e — ^ U^,, e 



-iEit 



aj 



M 



-iEit 



.i=i 



k=M+l 



2p 



;i-7) 



with Am^ = — ml, Ei = y p'^ + mf ^ p + mf/2p and t ^ L. Furthermore, using the 
unitarity relation 



M N 
j=l k=M+l 

we can rewrite the transition amplitude Au^^u^, in the form 



:i.8) 



Aua^l/0 — 6 



Ani^L , 
exp , 1 - 1 



N 



k=M+l 



1.9) 



From Eq.( |1.9| ), for the probability of —>■ vp transitions with /3 7^ a we obtain the following 
expression: 



_ 1 ( Am^L 

^Va^vp — 2 ^i^c-^yp (1 COS 2^ 



:i.io) 



where 



N 



A:=M+1 



:i.ii) 



is the amplitude of z/^ z/^ oscillations. 

The expression for the survival probability of can be obtained from Eq. ( [L.IUD and the 
conservation of the total probability: 



p = 1 - Vp 



1 



— 1 By I 1 — COS ■ 

" 2 °' " V 2p 



Arn^L 



:i.l2) 



where the oscillation amplitude By^.^^ is given by 



:i.i3) 
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Using the unitarity of the mixing matrix, from Eqs.( |l.ll| ) and (|1.13|) we obtain 



N 



N 



B =4 



\_k=M+l 



k=M+l 



Let us notice that from Eg. ( 1.11 ) we have 



A —A 



From Eqs. (|LTO|) and ( |LT5D it follows that 



P = P 



Due to CPT invariance we have (see, for example, Refs. 1 14, 36, 24 



P. 



= p- 

^up J- up- 



:i.i4) 



:i.i5) 



:i.i6) 



:i.i7) 



Thus, if only one squared-mass difference between the heaviest and lightest neutrinos is 
relevant for the oscillations of the terrestrial neutrinos, even in the case of non- conservation of 
CP in the lepton sector, the transition probabilities of neutrinos i^p) and antineutrinos 

{pet ^p) are equal: 



P = P- - 



:i.i8) 



Before closing this Section, let us remark that from the unitarity relation ( |1.8| ) it follows 
that the amplitude (|1.11| ) of Ua ^ z//3 oscillations can also be written as 



M 



a] 



:i.i9) 



In an analogous way, from Eqs. (|1.8| ) it follows that the oscillation amplitudes By^.,j^ given 
in Eq.( p..l4| ) can also be written as 



o=i / V i=i 



;i.20) 



Thus, the expressions for the oscillation amplitudes A^^-^^ and B^^.^y^ can be written with 
the sum over the indexes of the neutrinos of the second group (Eqs.( prTTD and ( |1.14| )) or 
with the sum over the indexes of the neutrinos of the first group (Eqs.(|lTT9D and ( |1.20| )). 



II. MIXING OF THREE MASSIVE NEUTRINOS 

Let us consider first the schemes with mixing of three massive neutrino fields. The case 
of a neutrino mass hierarchy 



mi <^ 777-2 <^ ms . 



(2.1) 
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was considered recently in Refs. p6| , p7| -|4^. In this case the oscillations of the terrestrial 
neutrinos are determined by three parameters: Am^ and the two mixing parameters |t/e3p 
and |?7^3p. In Ref. p8| , |39| it was shown that, if a hierarchy of couplings is realized in the 



lepton sector (i.e. |f/e3p and |t/^3p are small), the LSND result is not compatible with 
the results of all the other reactor and accelerator experiments on the search for neutrino 
oscillations. Instead, the result of the LSND experiment is compatible with the results of all 
the other experiments in the scheme with the neutrino mass hierarchy ( |2.1| ) if |f/e3p is small 
and |t/^3p is large (close to one). If massive neutrinos are Majorana particles, neutrinoless 
double-beta decay is allowed and for Am? > 5 eV^ this process can have a rate in the region 



of sensitivity of the next generation of experiments p6| , |39 |. 

Here we will consider another possible scheme with three massive neutrinos and one 
Am^ relevant for the oscillations of the terrestrial neutrinos. We will assume the following 
relation among the three neutrino masses: 



mi <^ 1712 — m-i 



(2.2) 

This scheme was recently discussed in Refs. pB|-pS|. In favor of such a scheme there are some 
cosmological arguments and some astrophysical arguments concerning the r-process 



production of heavy elements in the neutrino-heated ejecta of supernovae [35 



We will assume that the squared-mass difference A.m^2 



1712 is small and is 



relevant for the suppression of the flux of solar u^s on the earth. In this case, from the 
general formulas given in Section | it follows that the oscillations of the terrestrial neutrinos 
are characterized by the values of Am^ = m| — ml and two mixing parameters |t/eiP and 



|f/^i| . The probabilities of z/q, ^ {13 ^ a) and z/^ — > z/q, transitions are given by Eqs.( [L.lGD 



and ( |1.12| ) with the following oscillation amplitudes: 



4if//3inf/aii' , 

4|[/,ip(l-|f/„i|2) 



(2.3) 
(2.4) 



Let us consider first the constraints in the scheme under consideration following from 
the results of reactor and accelerator disappearance experiments. We will use the exclusion 
plots obtained in the z/g z/g Bugey reactor experiment and in the — >■ z/^ CDHS and 



CCFR84 accelerator experiments At fixed values of Am^, the allowed values of the 



oscillation amplitudes By^.^^^ and B^ are constrained by 



e,/i) 



(2.5) 



The values of 5°^.,^^ and B^^^-u^, ^e obtained from the corresponding exclusion curves. 
We will consider values of Am^ in the interval 



0.5 eV^ < Am^ < lO^eV^ 



(2.6) 



that covers the range where positive indications in favor of z/^ ^ Ue oscillations were reported 
by the LSND collaboration ||21[]. In this range of Am^ the quantities B^^.^^^ and B^^.^^ are 

smaller than 0.16 and 0.40, respectively. From Eq. (|2.4|) it follows that the parameters |?7qi|^ 
at fixed values of Am^ must satisfy one of the following inequalities: 
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|f/ai|'<a° (2.7) 



or 



with 



|f/ai|'>l-a° , (2.8) 



In the interval of Am^ given in Eq.( |2.6|) the values of and a° are smaller than 0.04 and 0.1, 
respectively. Thus, in the wide range of Am^ under consideration the mixing parameters 
\Ueif and I^T^il^ can either be small or large (close to one). 

We have assumed that Am32 is relevant for the suppression of the flux of solar neutrinos 
on the earth and that ^ Am^. In this case the survival probability of the solar u^s 

is given by 

P..^u. = {l-\U.ifYPi'M + \U.i\' , (2.10) 



(2 Si 

where P^e^u^ is the survival probability due to the mixing of with z/2 and z/3 (Eg. ( [2. 10] ) 



can be obtained with the method presented in Ref. [|^). From Eq.( |2.10|) and from the 
results of the Bugey disappearance experiment it follows that, in the case of a large |t/ei|^, 
we have Pu^-*v^ > 0.92 for all values of the neutrino energy. Such a large lower bound of 
the Ve survival probability is not compatible with the results of solar neutrino experiments 



1 2 1 2 

It is obvious that, due to the unitarity constraint \Uei\ + |f^ui| < 1, the mixing parame- 

I 1 2 I 1 2 

ters \Uei\ and \U^i\ cannot be both large. Thus, from the results of reactor and accelerator 
disappearance experiments and from the results of solar neutrino experiments it follows that 



1 2 1 2 

the values of the parameters \Uei\ and \Uui\ can lie in one of the following two regions: 



I. The region of small \Uei\^ and |f//^i|^. 
II. The region of small \Uei^ and large |f/^i|^. 

The fact that the parameter \Uei\' is small can have important implications for '^H (3- 
decay experiments and for the experiments on the search for neutrinoless double-beta decay 
((/?/?) oi^-decay). 

The electron spectrum in the decay ^He + e~ + z/g is given by 

^ = Cp, (g - T) F(E,) \Ue^? sj{Q-Tf-mj . (2.11) 

i 

Here Q is the energy release, Pe and Ef. are the electron momentum and energy, T = E^ — mf., 
E{Ee) is the Fermi function and C is a constant. Taking into account that \Ueif ^ 1 and 



7712 — 1TI3, from Eq. (|2.11| ) we have 



^ = Cp, E, {Q - T) F{E,) ^{Q-Tf-ml . (2.12) 
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This is the usual expression for the /5-decay spectrum from which information about the 
value of the "electron neutrino mass" m^^ is extracted. The upper limit for rriu^ obtained in 
the latest experiment is m^^ < 4.5 eV, which implies that 



Am^ < 20 eV^ 



(2.13) 



If z/2 and Us contribute to the dark matter in the universe, the masses m2 — are in the 
eV range. Thus, further improvements of the sensitivity of '^H /3-decay experiments could 
have important implications for the dark matter problem. 

If massive neutrinos are Majorana particles, neutrinoless double-beta decay of some even- 
even nuclei is allowed. The matrix element of {l3(3)ov decay is proportional to 



(m) = Ul 



nii 



(2.14) 



From the results of the experiments on the search for [(3j3)oy decay it follows that |(m) 



< 



1 — 2eV (see the review in Ref. |^). The expected sensitivity of experiments of the next 
generation is \ {m) \ ~ a few 10~^eV (see Ref. [pU[). 
In the model under consideration 

{m)^{Ul^ + Ul)m,. (2.15) 

We will consider the general case of non-conservation of CP in the lepton sector. The case 



of CP conservation was already discussed in Ref. p6| 
Let us write U^k as 



ek 



\Uek\ e 



If CP is conserved in the lepton sector we have |51 



Uek = ^efc Vk , 



(2.16) 



(2.17) 



where r/k = ±z is the CP parity of the Majorana neutrino with mass rrik- From Eqs.( p.l6D 
and (|2.17|) it follows that in the case of CP invariance we have 



±- 



TT 



(2.18) 



Furthermore, taking into account that |t/eiP 1, from the unitarity of the mixing matrix 
we have 



|f/e2|'+|t/e3| 

From Eqs. (gl5|) , (|2l6| ) and (^loD we obtain 



1 . 



|(m)| ^ m■i^J 1-4 |f/e2p If^esP siu^ $ 



(2.19) 



(2.20) 



with 
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$=103-021 . (2.21) 

If CP is conserved, we have $ = or 7r/2 for the same or opposite CP parities of P2 and z^s, 
respectively. 

Due to Eg. (|2.19|) we can write 

|?7e2|^cos6' and | f/es | ~ sin 6^ . (2.22) 

With the help of Eqs. (12^201 ) and (|]2|) we obtain 



\{m)\^ 



1 - sin^ 29 sin^ $ . (2.23) 



From Eq.(|]2|) it follows that 



l_sin2 20<M£<l. (2.24) 
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The boundary values of /rn^ correspond to CP conservation in the lepton sector: the 

upper (lower) bound corresponds to the case of equal (opposite) CP parities of z/2 and v^. 

Information about the value of the parameter sin^ 26 can be obtained from the results of 
solar neutrino experiments. If sin^ 2^ -C 1, which corresponds to the MSW solution with a 
small mixing angle, we have 

|(m)|^m3, (2.25) 

independently from the value of $ and from the conservation of CP. In the case of a large 
value of the parameter sin^ 26', which correspond to the MSW solution with a large mixing 
angle or to the vacuum oscillation solution, in the future it will be possible to obtain infor- 
mation about the violation of CP in the lepton sector if both the /9-decay and (/3/3)oi/ 
decay experiments will obtain positive results. In fact, from the measurement of in 
/3-decay experiments, of \{m) \ in {PP)ou decay experiments, and of sin^2^ in the solar neu- 
trino experiments, with the help of Eq. ( [^.23| ) it will be possible to determine the value of 
sin^ $ (if sin^ 29 is large). Let us emphasize that in the case of CP conservation the relative 
CP parities of z/2 and 1/3 can be determined: the CP parities are equal if sin^ <I> = and 
opposite if sin^ $ = 1. 

Now we will discuss the implications in the model under consideration of the results of 
the appearance neutrino oscillation experiments. We will use the exclusion plots, that were 



obtained in the BNL E776 experiment |^ searching for z/^ — > z/g transitions, in the FNAL 
E531 and CCFR95 experiments |]53| , |5^ searching for —>■ z/^ transitions and in the FNAL 



E531 experiment ||53| searching for z/g Vr transitions. We can use the results obtained 
in Ref. for the case of mixing of three neutrinos and the mass hierarchy ( p^.l| ) if we 

make the change 

\Ua.Z? -'Wal? . (2.26) 

In the following two subsections we will briefly present the main results for the two allowed 
regions of the values of the parameters |f/eiP and 



10 



A. The region | of small |J7eiP and \U^i 



2 



From Eq.( p.3|) it follows that the oscillation amplitude A^^-u^ is quadratic in the small 
quantities |?7eiP and whereas the oscillation amplitudes A^^-y^ and A^^.^^^ depend only 

linearly on these quantities. Thus, we expect that in the region | i/^ i/g oscillations are 
suppressed. 



From Eq.(2^), at fixed values of Am in the interval under consideration, we have 



A.,;.. <4a>;:, (2.27) 

with and a° given by Eq.( p.9| ). In Fig.|l] we have plotted the curve that represents this 
upper bound obtained from the results of the Bugey, CDHS and CCFR84 experiments 
(the curve passing through the circles). In Fig]^ we have also plotted the exclusion curves 
obtained in the BNL E776 |52| (dash-dotted fine) and KARMEN m (dash-dot-dotted line) 
experiments on the search of — >■ Ve transitions. The region allowed by the results of 
the LSND experiment is shown in Fig.^ as the shadowed region between the two solid lines. 
Taking into account that A^^.^y^ < By^.y^ (see Eq.( |1.13| )), we also plotted in Fig.|I|the exclusion 
curve for By^.y^ found in the Bugey experiment (dashed line). It can be seen from the figure 
that, for small values of Am^ (Am^ < 0.5 eV^), this bound on Ay^.y^ is stronger than the 
direct bound obtained by the BNL E776 and KARMEN experiments. It is also clear from 
the figure that this bound is not compatible with the result of the LSND experiment for 
Am^ < 0.2 eV^ 

It can be seen from Fig.|l] that, in the range of Am^ under consideration, with the 
exception of the region 10 eV^ < Am^ < 60 eV^, the limits on the oscillation amplitude 
A^^-^^ that can be obtained from the results of disappearance experiments are more stringent 
than the limits obtained in the direct experiments searching for z/^ —>■ z/g transitions. 

For Am^ > 4 eV^ we can obtain even stronger limits on the oscillation amplitude Ay^.^^ if 
we take into account the results of the FNAL E531 and CCFR95 experiments on the search 
for — > transitions. In fact, in the linear approximation over the small quantities |f/eiP 
and we have 

\U,A'^^. (2.28) 
From Eqs.( |2.3| ) and ( |2.28D we obtain the following upper bound: 

<«e<;.. , (2.29) 

where A^^.^^^ is the upper bound for the amplitude of ^ Vj- oscillations given by the 
exclusion curves obtained in the FNAL E531 and CCFR95 experiments. The bound (|2.29|) 



on the oscillation amplitude A^^.y^ obtained from the results of the Bugey, FNAL E531 and 
CCFR95 experiments is presented in Fig.|l] (the curve passing through the triangles). 

From Fig.|l] it can be seen that the results of reactor and accelerator disappearance 
experiments, together with the results of z/^ — appearance experiments allow us to 
exclude practically all the region of the parameters Am^ and Ay .y^ that is allowed by the 

1 2 1 2 

LSND experiment. Let us stress that this result is valid if the parameters \Uei\ and 
are both small. 
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In Fig.^ we have also shown the region in the Aiy^-iy^-Arn^ plane that will be explored 
when the projected sensitivity of the CHORUS and NOMAD |57| experiments, which 
are searching for z/^ — > transitions, is reached (the region limited by the line passing 
through the squares). 



B. The region |lj of small |?7eip and large |?7^i| 



From the discussion in Section [II A| it follows that, if the indications in favor of z/^ 



oscillations reported by the LSND collaboration are confirmed by future experiments, 
the parameters \Ueif and \U^if must satisfy the inequalities 

\Ueif < a° and \U^if > 1 - . (2.30) 

with the quantities 0° and a° given by Eg. ( p.9|) . 

Taking into account the unitarity of the mixing matrix and Eg. ( |2.30| ) we have 

\Ueif <l-\U,^f <al. (2.31) 



From Egs.( |2.30|) and ( p.31| ) we have 

\Ueif < Min [a°,aj] . (2.32) 



The inegualities ( p.30|) and the unitarity of the mixing matrix imply that also \Ut-i\ is small: 



|f/.i|'<l-|f/^i|'<aj . (2.33) 



Taking into account the inegualities ( |2.30D -( |2.33D , from Eg. ( |2.3[ ) it follows that the oscillation 



2 ^1 2 

amplitude A^^-^^ is guadratic in the small guantities \Uei\ and \Uri\ , whereas the oscillation 
amplitudes A^,^.^^ and A^^.y^ are linear in the same small guantities. For the oscillation 
amplitude A^^.^y,^ we have the following upper bound: 

A.,;., < 4 Min [a°,aj] aj . (2.34) 

In Fig.|^ we have plotted the upper bound for the oscillation amplitude A^^.^^ obtained from 
Eg. ( p.34| ) using the results of the Bugey, CDHS and CCFR84 experiments (the curve passing 



through the triangles). The solid line in Fig.|^ is the exclusion curve obtained in the FNAL 
E531 experiment. Taking into account the unitarity constraint A^^.^^ < By^-i,^, in Fig. ^ we 
have also plotted the exclusion curve obtained in the Bugey experiment. From Fig.|^ it can 
be seen that in all the considered range of Am^ the upper bound for Ay^.^^ given by Eg.( [^.34[ ) 
is very small, varying from about 10~^ to about 10~^ for Am^ > 0.3 eV^. 

Additional limits on the oscillation amplitude A^^.^^ can be obtained by taking into 
account the results of the experiments on the search for ^ z/g and z/^ ^ Vr oscillations. 
In the linear approximation in the small guantities |f/eiP, (1 — \U^ll\^) and If/riP, from 
Eg.(|2.3|) we have 

4 



f/eir^%^, (2.35) 



|f^^^|2^A^_ (2.36) 
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With the help of Eqs.( p.33D and ( |2.35|) we obtain the following upper bound for the oscillation 
amplitude A^^.^^: 

A.,,.^ < Al^.^^y^ . (2.37) 

In Fig.|^ we have plotted the corresponding boundary curve obtained from the results of the 
CDHS, CCFR84 and BNL E776 experiments (the curve passing through the squares). From 
this figure it can be seen that for Am^ > 1 eV^ the upper bound on the oscillation amplitude 
^ue;vr varies from 2 x 10^^ to 5 x 10^^ and is more stringent than that obtained from the 



results of disappearance experiments using Eq. (|2.34|) 
Finally, from Eqs.( p.35|) and (|2.36| ) we have 



A.,,.^ < '^""^ '^''"^ (2.38) 

The corresponding boundary curve obtained from the results of the BNL E776, FNAL E531 
and CCFR95 experiments is presented in Fig.|^ (the curve passing through the circles). 
From this figure it can be seen that for Am? > 4 eV^ the amplitude of ^ Ur oscillations 
is smaller than 3 x 10~^. 

Thus, if the parameter |t/eiP is small and the parameter is large, ^ J^r oscilla- 

tions are strongly suppressed. On the other hand, there are no constraints on the amplitudes 
of ^ Vf, and ^ v^. oscillations. As is well known, new experiments on the search for 
Vr transitions are under way at CERN (CHORUS J56[ and NOMAD fS^]). Experi- 
ments on the search for — >■ v^. transitions are continuing at Los Alamos (LSND [^) and 
at Rutherford Appleton Laboratory (KARMEN [Q). 

Let us conclude this Section with the following remarks. In the general case of mixing, 
the flavor neutrinos Vf.-, J^/x and Vr are not particles with a definite mass. Effective masses 
of the flavor neutrinos can be introduced in the framework of specific models. This is the 
case of the model under consideration. If the mixing parameters are in the region | i/g and 

are heavy and is the lightest neutrino: irij^^ ~ ms, m^^ ~ and niy^ ~ mi. If the 
mixing parameters are in the region |l| z/g and are heavy and is the lightest neutrino: 
m,^^ ~ m3, niij^ ~ ma and m,^^ ~ mi. Such inverted neutrino mass spectra are compatible 
with the r-process production of heavy elements in the neutrino-heated ejecta of supernovae 

m. 



III. MIXING OF FOUR MASSIVE NEUTRINOS 

All the existing indications in favor of neutrino masses and mixing (solar neutrinos, 
atmospheric neutrinos, dark matter, LSND) require at least three different scales of squared- 
mass differences. In this section we will consider short-baseline oscillations of the terrestrial 
neutrinos in a scheme with mixing of four massive neutrinos. We will consider a specific 



scheme ||29|-|34| with two light neutrinos z/i, z/2 and two neutrinos 1^3, z/4 with masses in the eV 
range. We will assume that the value of Am\.^ = ml — m\ is relevant for the suppression of 
the flux of solar i/g's and the value of Am2i = m\ — m\ could be relevant for the explanation 
of the atmospheric neutrino anomaly. Thus, we have 



13 



Am43 < Am^i < Am' 
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(3.1) 



with Am^ = m| — m^. 

The important difference between the model considered here and that considered in 
Section y is that in this model transitions of active neutrinos into sterile states are possible. 
The transition probabilities between different states are given by the general expressions 
( |l.l(j| ) and ( |1.12D with the following oscillation amplitudes: 



A 



4 



J=l,2 
4 6„(l-&a 



where 



i=i,2 



a] I 



(3.2) 
(3.3) 

(3.4) 



As in Section |I|, we start with the implications in the scheme under consideration of the 
results of the reactor and accelerator disappearance experiments. It is clear that we can use 
the results presented in Section We have 



or 



ha < a° (3.5) 

6a > 1 - < , (3.6) 

with a = e, /i. The quantities aP^ are determined by Eq.( |2.9| ). For each value of Am^ we 
determined the value of and from the exclusion plots obtained in the Bugey |H3 



CDHS @] and CCFR84 experiments. We will consider values of Am^ in the interval 
0.5 eV^ < Am^ < 10^ eV^, where the quantities and aj^ are smaller than 0.1. 

Thus, the values of the parameters h^. and can lie in one of the following four regions: 

I. The region of small and 6^. 

II. The region of small and large h^. 

III. The region of large hf. and small h^. 

IV. The region of large hf. and h^. 

If the neutrino masses satisfy the relation ( p.l|) and Am43 is relevant for the suppression 
of the solar flux, the probability of solar z/g's to survive is given by (see the Appendix ^ 



P. 



\ i=l,2 / i=l,2 



(3.7) 



where Pif^l^ is the survival probability due to the mixing of Ue with z/3 and z/4. If the 
parameter b^, is large, for the survival probability we have 
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Pu.^u. ^ E ■ (3-8) 

i=i,2 



Thus, if the parameters fog and he in the region |T| or in the region |^ the probabihty 
Pfe^ue practically does not depend on the neutrino energy. Moreover, taking into account 
that in the regions |T| and |^ we have \Uejf — 1, from Eq.( |3.8| ) we obtain 



The case of a constant Pu^-^u^ is disfavored by the solar neutrino data. In Ref. it was 
shown that the solar neutrino data cannot be explained with a constant Pp^^^^ if the fluxes 
of and ^Be neutrinos are allowed to vary in rather wide intervals around the standard 
model values. 

In both regions | and |T| we have \Uekf' — 1- Thus, information about the value 

fc=3,4 

of the masses ms ~ m4 can be obtained from the investigation of the end-point part of 
the /3-spectrum of decay (see Eq.( p.l2D ). If massive neutrinos are Majorana particles, 
neutrinoless double-beta decay is possible. It is clear that all the relations referring to this 
process derived in Section |I| in the framework of the scheme with three massive neutrinos 
(see Eqs.( p^.l3D -( p^.24D ) are also valid in the scheme under consideration. 



In the following we will discuss the implications of the results of neutrino oscillations 
appearance experiments if the values of the parameters and 6^ lie in the region | or in the 
region ||. 



A. The region ^ of small fee and fe^ 

In this region i/^ ^ z/g oscillations are suppressed. In fact, using the Cauchy-Schwarz 
inequality, from Eqs.( p.2|) and ( p.4D we obtain 

Au,,u. < Ah^he . (3.10) 

Taking into account that in this region < aP^ and 6e < 0,°, we obtain 

A.,;., <4aja°. (3.11) 

In Fig.^ we have plotted the curve that represent this upper bound obtained from the results 
of the Bugey, CDHS and CCFR84 experiments (the curve passing through the circles). 

Let us consider now the neutrino oscillation appearance experiments. From the unitarity 
relation ( p..l3| ) we have 

A.,;.. < 5.,;., , (3.12) 

V4.,;., < 5.,;., . (3.13) 

In the case of three massive neutrinos considered in Section H instead of these inequalities 
we had two approximate equalities, that allowed us to obtain additional strong limitations 
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for the amphtude of z/^ 



oscillations (see Eg. ( p.29|) ). In the model under consideration 



transitions of and into sterile states are possible and, instead of approximate equalities 
such as Eg. ( p .281 ), we have the inequalities 

^"^'^ , (3.14) 



be > 



A. 



(3.15) 



The inequality (|3.14]) can be useful if a positive result is found in the experiments on the 
search for z/^ ^ z/, oscillations (CHORUS [H, NOMAD fS^, COSMOS 11). 

An important difference between the model considered here and the one considered in 
Section |I| is that, in the case of mixing of four massive neutrinos, in the region | of small 
be and 6^ the result of the LSND experiment is compatible with the results of all the other 
neutrino oscillation experiments for 5eV^ < Am^ < 70 eV^. 



B. The region Jl| of small and large 6^ 

We will consider now the region fT|of small and large b^. From Eq.( |3.3D in the linear 
approximation in the small quantities be and 1 — 6^ we have 

be^^, (3.16) 

l-b^c^^^. (3.17) 

The essential difference of the model under consideration from the model with mixing of 
three neutrinos considered in the previous section is that, in the case of mixing of four 
massive neutrinos, Ue ^ z^r oscillations could be not suppressed in the region 0. This is due 
to the fact that in the model under consideration transitions of active neutrinos into sterile 
states are possible. In fact, using the Cauchy-Schwarz inequality, from Eq.(|3.2|) we have 



<46efer. (3.18) 

Furthermore, from the unitarity of the mixing matrix it follows that 

J2ba = 2. (3.19) 

a 

From this equation we cannot conclude that in the region of small be and 1 — 6^ the parameter 
br is small. Thus, the right-hand side of Eq. (|3.18|) could be not quadratic in small quantities 
(as it is in the case of mixing of three massive neutrinos) and Ai,^.^^ could be not suppressed. 

Let us notice that from the unitarity relation (|1.13| ) and from Eqs. (|3.16|) and (|3.17| ) we 
have the following lower bounds for be and I — b^ 



A 



be > ^ , (3.20) 
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C. Atmospheric neutrinos 



Up to now we have considered the constraints on the value of the mixing parameters be 
and due to the results of the solar neutrino experiments and of the reactor and accelerator 
neutrino oscillation experiments. If the squared-mass difference Am^i is relevant for the 
oscillation of atmospheric neutrinos (which is natural in this scheme), we can obtain an 
additional constraint on be and b^ from the results of the experiments on the detection of 
atmospheric neutrinos. 

In the model under consideration the survival probability of atmospheric is given by 
(see the Appendix |B]) 



p 



fl-Elf^-l') ^a..+ fEl^-l') ' (3-22) 

V fc=3,4 / \fc=3,4 / 

(I 2) 

where PuJ-!ua is the survival probability due to the mixing of with ui and z/2. In the 
region |, where both be and 6^ are small, we have 

P.^^u^ = bl Pl]:X + (1 - b,f ^ 1 , (3.23) 
Pu.^.. = blPllX + {^-hef^l. (3.24) 

Thus, if the parameters be and b^ are in the region |, the ratio of the observed neutrino- 
induced e-like and /z-like events in the atmospheric neutrino detectors should be equal to 



the expected ratio. This is in contradiction with the results of the Kamiokande |T8[, 1MB 



19| and Soudan experiments. On the other hand, if the parameters be and b^ are in the 



region ||, from Eq. ( p.22| ) we obtain 

Pu^,^u^, — Pul'-^u^ ) (3.25) 
Pu^-^u. ^ 1 , (3.26) 

and the results of the Kamiokande, 1MB and Soudan experiments can be explained with 
i^/j. ^ oscillations. 

Therefore, if the atmospheric neutrino anomaly found in the Kamiokande, 1MB and 
Soudan experiments is confirmed by future experiments, the values of the mixing parameters 
be and b^ can lie only in the region 



IV. CONCLUSIONS 

We have considered short-baseline neutrino oscillations, /3-decay and (/3/3)oy-decay in 
two specific schemes with mixing of three and four massive neutrino fields, respectively. 

Having in mind the experimental indications in favor of neutrino mass and mixing, we 
have assumed that in the case of three neutrinos the neutrino mass spectrum is composed 
of one very light neutrino with mass mi <C 1 eV and two neutrinos with masses m2 ~ in 
the eV range and a squared-mass difference m| — rn^ which is relevant for the suppression 
of the flux of solar i/g's. In this scheme, the oscillations of the terrestrial neutrinos are 
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determined by three parameters, Am^ = m\ — m\ and the two mixing parameters |f/eiP and 
From the analysis of the results of reactor and accelerator disappearance experiments 
and from the results of the solar neutrino experiments, it follows that only two regions of 
the values of the mixing parameters are allowed: |. |f/eiP and are both small; ||. 

|f/eiP is small and |f//iip is large (close to one). We have shown that in the region | the 
oscillations ^ z/g are suppressed and the result of the LSND experiment is not compatible 
with the results of all the other neutrino oscillation experiments. Instead, in the region |I| 
the oscillations z/g ^ z/^ are suppressed and we have obtained rather strong limits for the 
oscillation amplitude Ay^.^^^ from the results of neutrino oscillation experiments. If the 
masses m2 and are in the eV region, the effect of neutrino masses could be observed 
in the next generation of experiments on the measurement of the end-point part of the 
/3-decay spectrum and in the next generation of experiments on the search for {j3 (3) op-dec&y. 
We have shown that, if these experiment find a positive effect, a comparison of the results 
of the (/?/5)o,y-decay experiments and of the /3-decay experiments could allow to obtain 
direct information about the CP violation in the lepton sector. 

In the case of mixing of four massive neutrino fields we have assumed that the spectrum 
of neutrino masses is composed of two very light masses mi ~ m2 <C 1 eV with a value of 
m\ — m\ which could be relevant for the explanation of the atmospheric neutrino anomaly 
and two masses ~ m4 in the eV range with a value of m\ — m\ which is relevant for the 
suppression of the flux of solar z/g's. Since in this model transitions from active to sterile 
states are possible, the unitarity of the mixing matrix does not put as strong constraints on 
the oscillation channels as in the case of three neutrinos. The survival probabilities of z/g 
and z/^ are determined by the two parameters ha = If^aj T (with a = e, /i). From the 

results of neutrino oscillation disappearance experiments and the results of solar neutrino 
experiments it follows that only two regions of values of the parameters be and are allowed: 
|. The region of small be and 6^; |T[ The region of small be and large 6^. Oscillations z/^ ^ z/g 
are suppressed in the region |, but the LSND result is compatible with the results of other 
neutrino oscillation experiments if Am^ = m| — ml is in the range SeV^ < Am2 < 70 eV^ 
If the parameters be and b^ are in the region ||, z/g ^ u^. oscillations could be not suppressed 
(unlike the case of three neutrinos). If the atmospheric neutrino anomaly is confirmed by 
future experiments, the region | will be excluded. 

Both schemes considered here are in agreement with the cold and hot dark matter sce- 



nario with two practically degenerate neutrinos with masses in the eV range |32 



From the results of neutrino oscillation disappearance experiments and solar neutrino 
experiments it follows that in both schemes considered here the electron neutrino is "heavy". 
This means that both schemes are compatible with the constraints coming from the r-process 
nucleosynthesis in the neutrino-heated ejecta of supernovae ||35[| . 

In conclusion, let us remark that, if one of the neutrino mixing schemes considered here is 
realized in nature, i.e. there is no hierarchy of neutrino masses, the neutrino mass spectrum 
has a completely different character from the mass spectra of quarks and charged leptons. 
From our analysis of the results of neutrino oscillation experiments it follows that there is 
no hierarchy of couplings among generations in the lepton sector, i.e. the mixing matrix of 
neutrinos is completely different from that of quarks. 



18 



ACKNOWLEDGMENTS 

S.B. would like to acknowledge the kind hospitality of the Elementary Particle Sector of 
SISSA, where this work has been done. C.W.K. wishes to thank the Department of Physics, 
Korea Advanced Institute of Science and Technology, for the hospitality extended to him 
while this work was in progress. 



APPENDIX A: SURVIVAL PROBABILITY OF SOLAR NEUTRINOS 



In this appendix we will derive the formula ( p.7| ) for the survival probability Pu^^ue of 
solar neutrinos in the case of mixing of four massive neutrino fields considered in Section 

In general, the survival probability of solar u^s can be written as 



jk 



(Al) 



where f/(^^) is the mixing matrix in matter in the point of z/g production in the Sun and Ajk 
is the amplitude of z/™ — > transitions when the solar neutrinos travel from the central 
part of the Sun to the Earth, z/™ being the effective mass eigenstate neutrinos in matter in 
the point of production in the Sun (see, e.g., ||36|j2^ ). Equation (|A1|) is valid both for 



vacuum oscillations and for resonant MSW transitions of solar z/g's. In the case of vacuum 
oscillations f/^*^-* = t/, z/™ = uj and 



Ajk = Sjk exp 



—I 



Ej dx 



(A2) 



where Ej ~ p + m'j/2p, and Xi and Xf are the initial and final points of the neutrino 
trajectory. 

In the model under consideration it is assumed that Aml^ is relevant for the suppression 
of the flux of solar z/e's and Aml^ ^ -^^21 ^ Am^g. This implies that, if Am43 
range of the MSW solution of the solar neutrino problem, the transitions between z/^3^ and 
1^3(4) in the Sun can be strongly affected by solar matter effects. At the same time, the 
resonance densities associated with Amli ~ Am42 ~ Am^i ~ Am^g and with Am^i are 
much bigger than the density in the Sun and therefore the evolution of the z/i and 1^2 states 



in the Sun is adiabatic (see the discussions in Refs. pI|pD| ). As a consequence we have 



(M) 



U, 



for 



J = 1,2 



(A3) 



and Eq. (|A^) is valid for j', = 1, 2 and for j = 1,2 and k = 3, 4. Thus, the survival proba- 
bility averaged over the neutrino energy spectrum and the region of neutrino production in 
the Sun can be written as 



i=i,2 



fc,fc'=3,4 



(A4) 
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Using the unitarity of the mixing matrix, we can write Eg. as 

2 



where 



i=i,2 



ej I 



i=i,2 



I 



p(3,4) 



(A5) 



p(3,4) 



E ■ 

fc,fc'=3,4 



?7, 



(M) 



ek 



■■ -Akk' 



TT* 
'-'ek' 



fc"=3,4 



fc"=3,4 



In the case of vacuum oscillations Eq. ( |X6| ) can be written as 



p(3,4) 



1 - - sin^ 29 [l- cos 



2p 



(A6) 



which is the standard two-generation formula. Here the mixing angle 9 is defined by 



cos 9 



\U, 



e3| 



sin 9 



e4 



^fc=3,4 



(AT) 



(A8) 



\ \A:=3,4 



In the case of MSW resonant transitions Eq . can be written as 



P(3,4) ^ 1 , 



P34 ) COS 29m cos 29 



(A9) 



which has the form of the standard two-generation formula. Here P34 = I.4.34I is the 
probability of vf^ —>■ v^^ transitions, which is equal to the probability of i/™ transitions. 
The mixing angle 9 in vacuum is defined by Eq. (|A8|) and the mixing angle 9m in matter is 
defined by 



\U. 



cos 6* 



{M) 
e3 



\u. 



(M) 



M 



sin 9 



e4 



M 



(AlO) 



EK 



{M) 
ek 



\k=?,A 



Ek. 



(M) 
k 



=3,4 



APPENDIX B: SURVIVAL PROBABILITY OF ATMOSPHERIC NEUTRINOS 

In this appendix we will derive the formula ( |3.22| ) for the survival probability P^,^^^^ 
of atmospheric neutrinos in the case of mixing of four massive neutrino fields considered in 
Section |T|. 
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In general, the vacuum oscillation survival probability of z/^ is given by (see, for example, 



Refs. 14.36.2 



P^c^^yc ='^\Uaj\^\Uak\^ exp 



—I 



2p 



(Bl) 



with Am^j^. = m| — m^. 

In the model under consideration, for the atmospheric neutrinos we have 



exp —t 



2p 



1 . 



fB2) 



After averaging over the neutrino energy and the distance between the points of neutrino 
production and detection, we have 



exp —I 



2p 



~ for j = 1, 2 and A; = 3, 4 



(B3) 



From Eqs.(pTD-(p3D, for the survival probability of nUa we obtain 

2 



P. 



\A;=3,4 / j=l,2 



2 Am2i L 

a2 COS 

2p 



(B4) 



Using the unitarity relation (|1.8|) we have 



Elf. 

J=l,2 



aj I 



V fc=3,4 / 



2|t/ainf/, 



a2\ 



(B5) 



Combining Eqs.( p^ and (p5| ), we obtain the following expression for the survival probability, 

2 / \ 2 



P. 



\ fc=3,4 / \A;=3,4 / 



where 



U^2 


to 








2 \ 



1 — COS ■ 



Am2i L 



2p 



(B6) 



(B7) 



'aj I 



is the survival probability of due to its mixing with z/i and 1^2- In Section |I|we have used 
Eq.(il). 

Equation ( P7|) can be written in terms of a mixing angle 6^ defined by 



cos6'f, 



|f/al| 



sin 9a 



\U, 



a2\ 



(B8) 



, , E If. 

\ \i=i,2 
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Then the expression (p7| ) for the survival probabihty Pul'^ua takes the standard two- 
generation form, 

^£'i = l-^-n^2^^.(l-cos^f|^) (B9) 
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FIGURES 



FIG. 1. Exclusion regions in the A^^-^^-Atv? plane for small \Uei\'^ and |C//xiP in the model 
with mixing of three neutrinos discussed in Section The regions excluded by the BNL E776 and 
KARMEN —5- h'e appearance experiments are bounded by the dash-dotted and dash-dot-dotted 
curves, respectively. The dashed line represents the results of the Bugey experiment. The curve 
passing through the circles is obtained from the results of the Bugey, CDHS and CCFR84 exper- 



iments using Eg. ( 2. 27 ). The curve passing through the triangles is obtained from the results of 
the Bugey, FNAL E531 and CCFR95 experiments using Eq.( |2.29| ). The line passing through the 
squares bounds the region that will be explored by CHORUS and NOMAD. The region allowed 
by the LSND experiment is also shown as the shadowed region limited by the two solid curves. 

9 I 1 2 I 1 2 

FIG. 2. Exclusion regions in the A.i,^-i,^—Aui for small |Uei| and large |U^i| in the model 
with mixing of three neutrinos discussed in Section |ll[ The solid line represents the results of the 
FNAL E531 experiment. The dashed line represents the results of the Bugey experiment. The 
curve passing through the triangles is obtained from the results of the Bugey, CDHS and CCFR84 
experiments using Eq.( p.34 ). The curve passing through the squares is obtained from the results 



of the CDHS, CCFR84 and BNL E776 experiments using Eq.( ^) . The curve passing through 
the circles is obtained from the results of the BNL E776, FNAL E531 and CCFR95 experiments 
using Eq.(|]3|). 

FIG. 3. Exclusion regions in the A,y -j^^-Am? plane for small be and bn in the model with mixing 



of four neutrinos discussed in Section III. The regions excluded by the BNL E776 and KARMEN 
'^fj, t^e appearance experiments are bounded by the dash-dotted and dash-dot-dotted curves, 
respectively. The dashed line represents the results of the Bugey experiment. The curve passing 
through the circles is obtained from the results of the Bugey, CDHS and CCFR84 experiments 



using Eq.(3.11). 
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Figure |2| 
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